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A realistic estimate of the cusp effect in the rj — > 37r° decay is required for the forthcoming high 
precision experiments. The predictions for the size of this effect are given within the framework of 
nonrelativistic effective field theory. 
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I. INTRODUCTION 

The physical region in M n o„o invariant mass distribu- 
tion for r) — > 37T° decay extends below the charged two- 
pion threshold. It means that a cusp structure should 
be visible in this distribution around 2M 7r ± , in anal- 
ogy with the pronounced cusp in K + — > 7r + 7r 7r° decay, 
observed recently by NA48/2 collaboration [l[. In this 
paper, in particular, it has been shown that measuring 
charged kaon decays in the cusp region enables one to 
precisely determine S-wave tttt scattering lengths ao and 
(Z2, provided an accurate theoretical parameterization of 
the invariant mass distribution in terms of these scatter- 
ing lengths is known d !, i, G| (The strong impact of 
the unitarity cusp on tt°tt° scattering was already men- 
tioned in Ref. 0.). Moreover, the same logic applies 
to the neutral kaon decays into three pions, which have 
been studied in the recent experiment 0. The theoreti- 
cal framework for analysis of the neutral kaon decays is 
provided in Refs. [H, [HQ and the systematic inclusion of 
the electromagnetic effects both in charged and neutral 
kaon decays is considered in Ref. @ . We further mention 
that the general structure of the amplitude in the neutral 
kaon decays is similar to the rj — > 3ir° decay amplitude. 
For this reason, e.g., the two- loop representation of the 
amplitude in terms of the 7T7t effective-range expansion 
parameters, derived in Ref. [8], can be directly used to 
predict the cusp in the r\ — > 3ir° decay, which is studied 
in KLOE, Crystal Ball and WASA collaboration experi- 
ments Ha El Gi- 
lt should be pointed out that the two-loop formula 
for the kaon decay amplitudes, which was mentioned 
above, have been obtained in Refs. [H, Q within the 
non-relativistic effective field theory framework. This 
framework is ideally suited for parameterizing the final- 
state interactions in terms of the ttit scattering lengths 
(effective-range parameters, in general), whereas the ex- 
pansion of the amplitudes in Chiral Perturbation Theory 
(ChPT) is performed in powers of the quark masses and 
is less convenient for expressing the amplitude in the cusp 



region in terms of the observable quantities. (Note that, 
aside from the three-pion decays of charged and neutral 
kaons, the non-relativistic approach has been successfully 
applied recently to study of the K e 4 decays [l3|.) 

The aim of the present paper is to use the two-loop 
representation, derived in Ref. Q , to estimate the size of 
the cusp in the invariant mass distribution for 77 — ► 3ir° 
decays. This will finally allow one to judge, whether the 
forthcoming high-precision experiment will be able to see 
the cusp structure in the amplitude. Note that the cusp 
effect in the r\ — > 3ir° decay has been addressed already 
in various settings, e.g. in Refs. GHGH- 

In addition, we shall apply the same framework to 
study the experimental extraction of the slope parameter 
for the decay into three neutral pions. At present, the 
theory and experiment have not yet converged to a com- 
mon denominator for this parameter. ChPT at one loop 
in the isospin symmetry limit (l7j predicts a different sign 
for this parameter as compared to the experimentally 
measured one. At two loops, the sign of this quantity is 
no more fixed due to the large error bars coming from the 
unknown low-energy constants in ChPT [18j albeit the 
central value is still positive (The isospin-breaking correc- 
tions at one loop have been calculated in Refs. flq. llH |20| 
and are found to be small.). However, the predicted sign 
in Ref. [2l| where the calculations were done in the frame- 
work of unitarized ChPT, as well as the sign emerging in 
dispersive calculations [13, Hi, agree with the existing 
experimental data. We believe that in the forthcoming 
high-precision measurements of the slope parameter it 
will be very important to use as accurate a parameteri- 
zation of the decay amplitude, as possible. The param- 
eterization should be based on solid theoretical ground 
and, in particular, should take into account the cusp phe- 
nomenon which emerges at the physical values of the pion 
masses. 



II. THEORETICAL FRAMEWORK 

Below we mainly follow the notations from Ref. 
The tree-level amplitudes are expressed in terms of the 
kinetic energies Xi 
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where £Jj denote the pion energies in the eta rest frame. 
Up to the quadratic terms, 

Koo - Kv + K^Xl + Xt + Xl), 

M%% = U + L 1 X 3 + L 2 Xl + L 3 {X 1 -X 2 f , (2) 

where Li , Ki are the effective couplings in the non- 
relativistic Lagrangian that describe i] — > 37r decays at 
tree level. Note that we use the same notation for these 
couplings as in Ref. Q, where they denote the couplings 
describing the 3-pion decays of the neutral kaons. 

Assuming AI — 1 rule in the rj — > 3n vertex, the 
isospin symmetry relates the amplitudes for r\ — > 3tt° 
and r) — > 7T + 7r _ 7r° (we use Condon-Shortley phase con- 
vention) 

M 00 o(si,s 2 , s 3 ) = -M + - (si,s 2 , s 3 ) 

- M+-o(s2, S3, si) - M+-o(sa, sx, 8%) . (3) 

At tree level, this allows one to express the couplings Ki 
through Li 

K = -(3L + LxQ - L 3 Q 2 ) , 

K x = -(L 2 + 3L 3 ), (4) 

where Q = M v — 3M^o. 

In general, 77 — > 37r decay amplitudes are given in a 
form of a sum of the tree, one-loop, two-loop, . . . con- 
tributions Xqoo = + A^ loop + A4^ loops + • • • , 

and similarly for o- Thepertinent (rather lengthy) 

expressions are given in Ref. [8j . We do not display them 
here. It can be checked that these amplitudes in the 
isospin symmetry limit explicitly obey the constraints ([3|) 
at one- and two-loop level. 

We wish to stress that the representations given in 
Refs. H [s] should be understood as a parameterization 
which should be fit to the data. In other words, the con- 
stants Li, do, a 2 , • • ■ are considered as free parameters to 
be fixed from the fit. In this paper, we however make an 
attempt to predict the size of the cusp - fitting first the 
tree-level amplitude in order to determine Li and then 
using one- and two-loop representation to produce the 
cusp in the synthetic data. In doing this, we have fixed 
ao, a,2 to their theoretical values [24[ and neglected isospin 
breaking in the derivative 4-pion couplings, as well as the 
shape parameter and the P-waves. 

The matching of Li is done to: 

1) The tree-level amplitude 77 — * 7r + 7r~7r in ChPT. 
According to Eq. the overall normalization of 
the amplitude does not play a role, only the slopes 
matter. The result is given by 

L = (4M2 -3(A^-A^o)) 2 /(Af2-M2 ), 
Lx = 6M V /(M*-M% ), 

L 2 = L 3 = 0. (5) 
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FIG. 1: Invariant mass distribution dT/dM n o^o divided by 
the phase space, calculated at two loops: 1) Matching to 
ChPT at tree levelfdashed line); 2) Matching to the KLOE 
parameterization [251 ] (solid line). 

2) The experimental amplitude extracted by KLOE 
collaboration [25| . 

In order to carry out the matching to the KLOE data, 
it is useful to introduce Dalitz variables for 7/ — > 7r + 7r~7r° 
decay 

x = V3{Xx-X 2 )/Q, y = 3X 3 /Q-l. (6) 

For the decay r\ — > 3ir° one defines the variable 

z = x 2 + y 2 . (7) 

The phenomenological parameterization of the amplitude 
is given by 

M+- a = A c (l + ay + f3y 2 + -/x 2 ) , (8) 

with a, (3, 7 being complex quantities. The matching of 
Eq. © to the real part of Eq. © yields 

L = A c (l - Re a + Re/3) , 

Lx = 3A c (Rea- 2Re/3)/Q, 

L 2 = 9A c Re(3/Q 2 , 

L 3 = 3A c Rc 7 /Q 2 . (9) 

The right-hand side in Eq. Q is fixed by using Eq. (6.4) 
and table I of Ref. [25] . Isospin-breaking corrections in 
Eqs. ([8]) and ([9]) are consistently neglected. 

We would like to mention that the systematic way of 
fixing the parameters of the effective non-relativistic La- 
grangian consists in performing a simultaneous fit of the 
non-relativistic representation to both charged and neu- 
tral invariant mass distributions. The results, which are 
contained in the present paper, should be considered only 
as a rough theoretical estimate of the expected size of the 
cusp effect in the 77 — > 37r° decay. 



3 




0.28 0.30 0.32 0.34 0.36 0.38 0.40 

Mn°7i° [GeV/c 2 ] 

FIG. 2: Invariant mass distribution dT/dM n o^o divided by 
the phase space, with the couplings matched to the KLOE 
parameterization !25j: 1) Tree level (dotted line); 2) One loop 
(dashed line); 3) Two loops (solid line); 4) Two loops, assum- 
ing M n ± = M n o (dot-dashed line). 
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FIG. 3: Differential decay rate dT/dz divided by the phase 
space at two loop: I) Matching to ChPT at tree level (dashed 
line); 2) Matching to the KLOE parameterization [25j . 



III. RESULTS 

In Fig. [T] we display the calculated invariant mass dis- 
tribution for rj — > 3ir° decay, divided by the phase space. 
The decay amplitude is normalized in the center of the 
Dalitz plot 



M 2 

7) 



|A4ooo(so,s ,so)| =1, s = + M. 



(10) 



We display the result for Li, Ki matched to the tree-level 
result of ChPT, or to the KLOE amplitude [25]. The 
resulting cusp in both cases amounts roughly up to a 2% 
effect. We would like to mention that the sign of the cusp 
effect is fixed by the isospin symmetry, see Eqs. ([3]) and 
(|3]) and is thus a robust theoretical prediction. 

In order to check the convergence of the method, in 
Fig. (J2) we show the invariant mass distribution calcu- 
lated at tree level, one and two loops, with the couplings 
Li matched to the KLOE amplitude. It is seen that the 
shape of the cusp does not change much from one- to two- 
loop calculations, indicating at a rather robust prediction 
for a size of this effect. 

Figure [3] contains our prediction for the differential de- 
cay rate in the variable z - again with Li , Ki matched 
either to the tree-level result of ChPT, or to the KLOE 
amplitude. As expected, the slope parameter in the for- 
mer case has the opposite sign as compared to the ex- 
perimentally observed. Apart from a small dip around 
z ~ 0.75, corresponding to the cusp, the differential de- 
cay rate is seen to be fairly linear in z. 

The convergence of the loop expansion for the differen- 
tial decay rate in the variable z and the comparison with 
the Crystal Ball data [l(J is shown in Fig. 2J As seen, 
the data are described quite well. 



IV. CONCLUSIONS 
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FIG. 4: Calculated dV/dz divided by the phase space and 
comparison with Crystal Ball data [To(] . Couplings matched 
to the KLOE parameterization [25[. 1) Tree level (dotted 
line); 2) One loop (dashed line); 3) Two loops (solid line). 



Using the two-loop parameterization of the 77 — ^ 37r 
decay amplitudes [1| , we have shown that the size of the 
cusp effect in the invariant mass distribution for i] —>■ 3ir° 
process amounts up to around 2%. Despite such tiny ef- 
fect, one may expect that forthcoming high-precision ex- 
periments at Crystal Ball, KLOE and WASA-at-COSY 
with about 10 7 events in the Dalitz plot will be able to 
observe it. It is however unlikely that one could deter- 
mine 7T7r scattering lengths at a reasonable accuracy from 
these experiments. 

Moreover, the cusp effect modifies the differential de- 
cay rate for the r\ — > 37r° decay in the variable z, produc- 
ing a dip around the value z ~ 0.75. We expect that, in 
order to carry out an accurate analysis of the Dalitz plot 
distributions, this effect should be taken into account. 

Finally, we wish to mention that for the cusp-like struc- 
ture, which has been seen recently by the Crystal Ball 
collaboration experiment at MAMI-C (26j, the sign of 
the effect is claimed to be different from the theoretical 
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prediction. To resolve this contradiction, experimental 
study of the r\ — > 3ir° decay with a better statistics would 
be desirable. 
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